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SUMMARY 

The  limits  of  applicability  of  the  quasi-steady  state 
assumption  on  which  the  modified  shell  trajectory  model  is 
based  have  been  investigated.  Conclusions  have  been  drawn 
about  the  applicability  of  the  modified  point  mass  model 
to  fire  prediction  for  very  high  gun  elevations. 
Implications  for  calibration  of  shell  fire  control  models 
using  the  Mach  number  fitting  process  are  particularly 
interesting. 
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,1.  INTRODUCTION 


Modern  computerised  methods  for  ^reducing  ordnance  firing  tables  are  based  on 
the  modified  point  mass  trajectory  model  for  a  spin  stablised 
projectile(ref . 1) .  The  model  includes  only  the  effects  of  the  steady-state 
gravity- induced  trim  or  yaw  of  repose.  In  this  way  significant  contributions 
from  the  angular  motion  of  the  projectile  are  included  in  the  trajectory 
computation,  but  the  penalties  arising  from  the  very  small  integration  step 
size  necessary  for  accurate  representation  of  the  transient  angular  motion  are 
avoided.  The  modified  point  mass  model'  is  valid  for  a  dynamically  stable 
missile  and  slowly  varying  yaw  of  repose.  ^4lurphy(ref . 2)  has  shown  that  for  a 
sufficiently  high  angle  of  fire  the  quasi-steady  state  assumption  breaks  down 
near  apogee,  where  conditions  change  most  rapidly.  When  this  occurs 
precession  or  the  slow  mode  of  the  epicyclic  motion  is  excited  and  transient 
motion  persists  on  the  downleg.  As  a  result  trajectory  predictions  with  the 
modified  point  mass  model  show  significant  range  errors  for  some  high  angle 
trajectories . 


nurphy  suggested  in  reference  2  that  a  further  modification  could  be 
incorporated  in  the  modified  point  mass  model  to  correct  this  problem.  The 
idea  has  not  been  followed  up,  principally  because  the  gun  elevations  which 
produce  such  behaviour  are  not  used  a  great  deal  in  practice.  However,  the 
break  down  of  the  quasi-steady  assumption  does  have  implications  for  the 
calibration  of  the  modified  point  mass  trajectory  model  for  shell  fire  control 
purposes.  The  modern  Mach  number  fitting  technique  for  calibration  of  shell 
fire  control  models (ref . 3) ,  using  data  from  Range  and  Accuracy  trials,  relies 
on  the  overall  accuracy  of  the  trajectory  model,  Severe  distortions  may  arise 
in  the  calibration  process,  and  in  using  the  tfiodel  for  prediction  purposes  if 
there  are  significant  errors  in  trajectory  predictions  for  high  angle  fire. 


ftiis  brief  report  summarises  the  results  in  reference  2  and  looks  at  the 
magnitude  of  the  transient  motion  induced  by  the  rapidly  changing  conditions 
near  apogee  and  the  relation  between  that  motion  and  gun  elevation.  Results 
are  given  for  a  variety  of  shells  of  interest  to  the  Australian  Army.  Finally 
some  conclusions  are  drawn  about  the  relevance  of  the  effect  to  practical 
trajectory  calculation,  with  particular  emphasis  on  model  calibration  and 
Range  and  Accuracy  trials . 


2 .  THEORY 


Two  right  handed  orthogonal  sets  of  axes  are  used  in  describing  projectile 
motion.  The  earth  fixed  axes  have  the  x_-axis  downrange,  z  -axis  vertically 


e 


downwards  and  ye-axis  to  the  right, 


The  aeroballistic  axes  have  the  x-axis 


along  the  axis  of  symmetry  of  the  projectile,  the  z-axis  in  the  plane  of  the 
trajectory  pointing  downwards  and  the  y-axis  horizontal  to  the  right.  The 
complex  angle  of  attack,  defined  in  the  aeroballistic  axes  system  is 


A  A  A 

K  =  (v  +  iw)/V , 


wher  v,  w  are  y  and  z  velocity  components,  so  that  the  real  part  of  5  is  in  a 
horizontal  plane  and  the  imaginary  part  is  in  a  vertical  plane.  A  detailed 
derivation  of  the  gravity- induced  yaw  of  repose  for  a  spin  stabilised 
projectile  is  given  in  reference  1.  A  very  good  approximation  for  the  yaw  of 
repose(ref .2)  is 


Pi„  sz/Qsdvca 


r*j| 


a 


-V,| 


Oi 


CD 
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where  p  =  roll  rate, 

I  =  axial  moment  of  inertia, 

g  =  gravity  component  along  z-axis, 
z 

Q  =  dynamic  pressure, 

S  =  projectile  cross-sectional  area, 
d  =  projectile  calibre, 

V  =  relative  air  velocity  of  projectile,  and 

C  =  overturning  moment  coefficient, 
a 

If  the  yaw  of  repose  is  rapidly  varying  near  apogee  the  total  yaw  can  be 
represented  by 

5  =  -6  +  6  K_  (5  ,  02)  e1*2  (2) 

g  g  ga  G  ga  ' 

where  6  is  the  value  of  6  at  apogee,  and 

ga  g 

02  is  the  phase  angle  of  the  slow  mode  of  the  epicyclic  motion  and 

damping  of  the  slow  mode  is  ignored  for  simplicity.  The  phase  angle 

can  be  derived  by  integrating  the  equation 

02  =  (plx/ly)  [1  -  Cl-Sg"1)*]  (3) 

where  s  -  (pi  )2/4QSd  I  C  is  the  gyroscopic  stability  factor.  The 

g  x  y  m0 

amplitude  function  Kr(6  ,02)  is  defined  by  the  relations 


KG  =  °» 


02  S  "2tr,  upleg 


KG  -  Sga 


f(0.p)  d02,  -2ir  <  02  <  2ir,  near  apogee 


where 


kg  "  1  K2q  (6ga)  **  *  2ir»  down  leg 


r 


f(0_)  sin#2 1  d  ♦«* 


and  the  function  f(8^,)  is  defined  by  the  relations 

f(0T)  =  4  cos7  0T  sin  0T  (5) 

36ga  =  -tan  0T  (3  +  tan4  0T)  (6) 

since  0^  is  the  angle  between  the  trajectory  and  its  projection  onto  the 
horizontal  plane.  The  zero  for  the  phase  angle  $2  is  chosen  to  coincide  with 
apogee,  when  8^=0,  so  that  $2  <  0,  relates  to  the  upleg  of  the  trajectory  and 
>  0  to  the  downleg.  Then  the  total  yaw  on  the  downleg  of  the  trajectory  is 


<G  -  -Sg  +  1  6ga  K2G  <  V  '  <7> 

The  significant  effect  of  the  increased  yaw  on  the  trajectory  is  to  increase 
the  yaw  drag.  The  overall  drift  will  not  be  substantially  changed,  because 
only  averaged  results  affect  the  impact  point  and  the  averaged  effects  of  the 
change  in  lift  will  be  small.  The  total  drag  coefficient  is  defined  as 


CD  “  CD  +  CD  ,  V  •  <»> 

o  a 

On  the  upleg  of  the  trajectory  Kq=0,  so  that  the  modified  point  mass  model  is 

accurate.  Near  apogee,  the  projectile  velocity  is  low  and  small  changes  in 
drag  do  not  affect  the  trajectory  significantly.  However,  on  the  downleg,  the 
effects  can  be  significant.  The  object  of  this  paper  is  to  establish  the 
conditions  under  which  the  effect  may  be  significant  and  to  look  at  methods 
for  dealing  with  the  resulting  problems  in  using  the  modified  point  mass  model 
for  fire  prediction. 


3.  NUMERICAL  RESULTS 

The  first  step  in  evaluating  the  significance  of  the  effect  is  to  determine 
the  amplitude  excited  in  the  slow  mode  of  the  motion.  Figure  1  shows  how  the 
function  X^q  (6gft)  varies  with  yaw  of  repose  at  apogee  and  the  resultant 

amplitude  of  the  epicyclic  motion  which  is  excited.  The  initial  amplitude  of 
the  epicyclic  motion  has  a  negative  maximum  of  about  0.75°  for  a  yaw  of  repose 
at  apogee  of  about  6°.  The  amplitude  then  changes  sign  and  the  magnitude  does 
not  reach  1®  until  yaw  of  repose  is  almost  10°,  after  that  it  increases 
rapidly.  Figure  2  shows  a  typical  example  of  the  motion  excited.  It  depicts 
the  incidence  history  for  an  Ml  shell  launched  at  an  elevation  of  70°.  The 
slow  mode  motion  on  the  downleg  is  typical  of  this  sort  of  trajectory. 

In  order  to  evaluate  the  consequences  of  these  results  let  us  look  at  how  yaw 
of  repose  varies  with  gun  elevation  for  a  variety  of  shells.  Figure  3  shows 
yaw  of  repose  for  two  105  mm  projectiles  and  two  155  mm  projectiles.  Results 
are  plotted  for  muzzle  velocities  equivalent  to  maximum  charge  for  each 
projectile.  The  yaw  of  repose  for  which  the  induced  epicyclic  motion  has  zero 
amplitude,  8.5°,  is  achieved  for  elevations  between  63®  and  67®.  For  a  gun 
elevation  of  70°,  yaw  of  repose  at  apogee  varies  from  14.5°  for  the  Ml  shell 
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to  over  30°  for  the  L15.  Therefore  some  shells  will  exhibit  significant 
excitation  of  the  slow  mode  of  their  epicyclic  motion  for  an  elevation  of  65°, 
and  all  shells  can  be  expected  to  show  the  effect  for  an  elevation  of  70°. 
Figure  4  shows  the  range  error  of  the  modified  point  mass  model  for  the  Ml 
shell ,  which  has  the  slowest  development  of  the  phenomenon.  The  yaw  of  repose 
of  the  Ml  shell  at  apogee  is  8.5°  for  a  gun  elevation  of  66.6°.  Above  this 
elevation  it  is  clear  from  figure  4  that  the  range  error  in  modified  point 
mass  model  predictions  increases  very  rapidly.  We  can  therefore  regard  the 
elevation  at  which  yaw  of  repose  at  apogee  is  8.5°  as  a  critical  value  for  the 
modified  point  mass  model,  since  the  accuracy  of  model  predictions 
deteriorates  rapidly,  from  that  point  onwards.  Table  1  shows  the  elevations  at 
which  this  critical  value  is  reached  for  each  of  the  shells  represented  in 
figure  2.  Firing  tables  generally  give  data  for  elevations  up  to  70°  and 
consequently  model  calibration  uses  data  from  Range  and  Accuracy  firings  which 
contains  significant  components  for  elevations  around  70°.  Therefore  this 
effect  will  produce  significant  distortion  of  the  calibration  and  fire 
prediction  process  unless  handled  carefully. 


I 


H 


El 


4.  CONCLUSIONS 

It  is  clear  that  for  most  shells  the  flaw  in  the  modified  point  mass  model 
which  has  been  discussed  will  not  have  any  significant  effect  for  elevations 
below  65°.  Since  the  need  to  use  artillery  at  higher  elevations  that  this  is 
extremely  rare,  it  appears  at  first  glance  that  the  effect  can  be  quite  safely 
ignored.  However,  problems  do  arise  in  the  production  of  firing  tables.  In 
general,  firing  tables  are  produced  for  elevations  up  to  70°,  and  therefore 
the  Range  and  Accuracy  trials  on  which  the  tables  are  based  often  contain  up 
to  20%  of  data  points  from  firings  at  elevations  between  65°  and  70°,  and 
sometimes  even  more.  Because  of  the  rapidly  deteriorating  performance  of  the 
modified  point  mass  trajectory  model  at  these  elevations,  the  calibration 
process  has  to  cope  with  significant  errors  from  this  part  of  the  data  which 
will  produce  corresponding  distortions  of  the  fitting  process  at  low  to 
moderate  elevations.  The  effects  will  be  particularly  severe  for  the  yaw  drag 
calibration  factor,  QFAC,  in  the  Mach  number  fitting  method(ref . 3)  because  of 
the  global  nature  of  the  process. 

The  alternative  suggested  in  reference  3  of  a  revision  of  the  modified  point 
mass  scale  would  probably  improve  the  trajectory  representation  sufficiently 
to  avoid  distortion  and  bias  effects.  However,  since  the  firing  tables  are 
used  so  rarely  at  these  high  elevations  it  is  possible  to  restrict  the 
ballistic  data  gathering  trials  and  the  follow  up  .calibration  process  to 
elevations  for  which  the  yaw  of  repose  at  apogee  is  below  the  critical  value 
of  8.5°.  This  would  form  an  additional  constraint  on  the  design  method  for 
Range  and  Accuracy  Trials,  described  in  reference  4. 


TABLE  1.  CRITCAL  GUN  ELEVATIONS  FOR  6  =8.5° 

ga 


Projectile 

Calibre 

(mm) 

MV 

(m/s) 

Elevation 

(°) 

Ml 

105 

464.8 

66.6 

L31 

105 

708.0 

65.2 

L15 

155 

826.2 

63.5 

M483A1 

155 

643.0 

64.7 
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